Microtubule length distributions in the presence of protein-induced severing 
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Microtubules are highly regulated dynamic elements of the cytoskeleton of eukaryotic cells. One 
of the regulation mechanisms observed in living cells is the severing by the proteins katanin and 
spastin. We introduce a model for the dynamics of microtubules in the presence of randomly 
occurring severing events. Under the biologically motivated assumption that the newly created plus 
end undergoes a catastrophe, we investigate the steady state length distribution. We show that the 
presence of severing does not affect the number of microtubules, regardless of the distribution of 
severing events. In the special case in which the microtubules cannot recover from the depolymerizing 
state (no rescue events) we derive an analytical expression for the length distribution. In the general 
case we transform the problem into a single ODE that is solved numerically. 

PACS numbers: 87.16.Ka, 87.16.ad, 87.10.Ed 



I. INTRODUCTION 

Microtubules are filamentous protein aggregates that 
appear in all eukaryotic cells. They have an inherent po- 
larity that results in different dynamics at their two ends. 
The so-called plus end is highly dynamic, alternating be- 
tween prolonged periods of polymerization (growth) and 
dcpolymcrization (shrinkage) [1]. On the other end of 
the filament, the minus end often remains connected to 
the locus of nucleation [2] or is found to exhibit rela- 
tively steady depolymerization [3]. The combination of 
slow depolymerization at the minus end and prolonged 
growth at the plus end leads to a phenomenon known 
as treadmilling, whereby the individual tubulin dimcrs 
appear to move from the plus to the minus end [4]. As 
the stiffest of the cytoskeletal filaments, microtubules are 
widely used in intracellular transport and for structural 
support. Therefore, the dynamic properties of the micro- 
tubules are tightly regulated by the cell, often through 
the use of microtubule-associated proteins (MAPs) [5] . 

In this work we investigate one particular type of MAP 
that causes the severing of microtubules. It was noted by 
Vale [6] that otherwise stable microtubules could be sev- 
ered in mitotic extracts of Xenopus eggs. This activity 
was traced back to a protein that is able to use ATP 
hydrolysis to sever microtubules. The protein was iden- 
tified only later and given the name katanin after the 
katana, the Japanese Samurai sword [7]. Katanin is a 
hcterodimer, consisting of the p60 and p80 subunits: the 
p80 subunit is thought to be responsible for the target- 
ing of the protein, whereas p60 is involved in the actual 
severing as part of a hexameric ring [8]. The hexameric 
form of katanin appears to remove individual dimcrs from 
the microtubule lattice, thereby compromising the struc- 
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tural integrity of the microtubule. It is not currently clear 
whether katanin acts uniformly along the microtubule, or 
whether it is attracted by pre-existing lattice defects [9]. 

Katanin homologs have since been discovered across 
the animal and plant kingdoms [2]. Another severing 
protein by the name of spastin has also been identified. 
Like katanin, it also assembles in a hexameric ring, sug- 
gesting a severing mechanism similar to that of katanin 
[10]. Both katanin and spastin are capable of severing mi- 
crotubules at seemingly random locations, and are used 
by the cell for the regulation of the cytoskeleton [11], for 
example in the mitotic and meiotic spindles [12] or the 
formation of the transverse cortical array in plant cells 
[13, 14]. Generally, the activity of severing proteins leads 
to a decrease in the average microtubule length, but an 
increase in their number [11]. Surprisingly, this increase 
in number can sometimes more than offset the loss of 
microtubule length due to the average length decrease 
[15]. 

In this work, we investigate theoretically how the oc- 
currence of microtubule severing at random positions af- 
fects the length distribution of microtubules. Previous 
studies have assessed the effect of severing on actin fil- 
aments [16, 17]. However, whereas actin has a single 
growth mode that is described well by constant polymer- 
ization and depolymerization rates, microtubules show 
richer dynamics because the plus end switches between 
growing (polymerizing) and shrinking (depolymerizing) 
states. 



II. MODEL 

We base our model on the basic dynamic instability 
model that was introduced by Dogterom and Leibler [18]. 
In this model, microtubules exist in either the growing 
or the shrinking state, and switch between these states 
with a 'catastrophe' rate r c (growth — > shrinkage) and 
a 'rescue' rate r r (shrinkage — > growth). In the growing 
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state, the microtubule extends with an effective speed 
v + , and in the shrinking state it retreats with an effec- 
tive speed v~ . New microtubules are nucleated with a 
steady nucleation rate r n . This set of constraints gives 
rise to a steady state distribution of microtubule lengths, 
provided that v + /r c < v~ /ry. 

We model microtubule severing by a constant severing 
rate per unit length. This is a valid approximation if 
we assume that microtubule severing process occurs on 
a time scale that is much shorter than the time in which 
a microtubule grows significantly. By taking a constant 
severing rate, we also implicitly assume that severing is 
limited by the availability of microtubules, i.e. the sever- 
ing protein is available in abundance. However, because 
we focus on steady-state results, where the total amount 
of microtubules is constant, the possible invalidity of this 
assumption will not qualitatively affect the results. 

Following the approach in [18], we construct a set of 
evolution equations for the length distributions of grow- 
ing and shrinking microtubules. Denoting the grow- 
ing and shrinking microtubule length distributions by 
m + (l,t) and m~(l,t), respectively, the equations can be 
written as 
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where the derivatives with respect to / reflect the trans- 
lation of the distributions due to growth and shrinkage 
of microtubules. The severing contribution $ 
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be constructed below. These equations are supplemented 
by the boundary condition 



■>+' 



(2) 



specifying the nucleation of new microtubules with rate 
r n , and by the physically motivated constraint that both 
distributions tend to zero for large lengths (there are no 
infinitely long microtubules). 

In the following we derive the form of the severing 
terms in the set of equations (1). The inclusion of sev- 
ering events leads to two types of contributions to the 
evolution of the length distributions: the disappearance 
flux </w(M) °f microtubules of a certain length I and 
the appearance fluxes <pi n (l',t) and </>;„(/", t) of two new 
microtubules with a total length I = I' + I" . In addition, 
the newly created microtubules can be created from mi- 
crotubules that were initially growing (+) or shrinking 
(-). Symbolically, we write 

t (M) + 0+ + (M) + 0+_(M) (3) 

^outM + C,+ (M) + fe,_0,i), (4) 

where (f>^ nt (l,t) stands for the removal of microtubules 
in state a £ { — , +}, and 4>f nT (l, t) represents the appear- 
ance of new microtubules of length I in state a as a result 
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Table I: Overview of the parameters in natural dimensions 
and the dimensionless parameters 



of the severing of microtubules in state r. The contribu- 
tions will be discussed individually below. 

The process of severing is controlled by the severing 
rate r s that is given as a rate per unit of length. Thus we 
find that the fluxes of disappearing microtubules 4>^ ut (l, t) 
(growing) and </>~ ut (/,i) (shrinking) arc given by 



0out(M) 

<feut(M) 



r s lm + {l,t) 
r s lm~ (I, t). 
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In order to specify the influx terms 



5 in ± it is necessary 
to specify the process of severing in more detail. We will 
assume that the action of the severing protein is local, 
thus having no effect on the remote plus and minus ends 
of the microtubule it severs. This implies that the plus 
end fragment of a severed microtubule remains in the 
same state. However, we must make an explicit assump- 
tion regarding the state of the newly created plus end. In 
line with biological observations [19], we assume that this 
plus end always starts out in the shrinking state. A sever- 
ing event thus shortens an existing microtubule without 
affecting its growth state and creates an additional mi- 
crotubule that is in the shrinking state. If necessary, the 
model can easily be extended to handle (a fraction of) 
severing-created plus ends that start out in the growing 
state. 

Let us define the uniform probability density of select- 
ing a severing location I on a microtubule of length 
as p(l\l'), with the value l/l' for I £ [0,1% and oth- 
erwise. We then derive the influx t) of growing 
microtubules of length I that results from the severing of 
growing microtubules. 
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where 6(x) is the Heavisidc step function. In a similar 
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way, and taking account the fact that the minus end frag- 
ment is always in a shrinking state, we derive 



<t>i ut (i',t)p(i\i')<u' 



r s J m+(l')dl' 



(8) 
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2r s J m-{l')dl' 



(9) 



And, because the severing of shrinking microtubules can- 
not produce growing microtubules, we have (f)f n _{l,t) = 
0. 

Note that we do not need to keep track of the correla- 
tions between the lengths of the individual microtubules 
that are created from a single cutting event, because we 
are looking only at ensemble-averaged length distribu- 
tions. The contributions derived above are similar in 
form to those introduced by Edelstein-Keshet and Ermen- 
trout [16] as the continuous limit of a discrete monomer 
addition and severing model for actin filaments. How- 
ever, we use distinct asymmetric terms for growing and 
shrinking microtubules. 

The model as described does not explicitly take into 
account the often observed depolymerization of micro- 
tubules at their minus ends that leads to treadmilling. 
However, this case is easily addressed through a renormal- 
ization of the growth and shrinkage speeds. Denoting the 
shrinkage speed at the minus end by v tm , this is achieved 



by the substitutions it 



and v 



leaving the results qualitatively unchanged by the incor- 
poration of treadmilling. 



A. Steady state equations 

Inserting the severing terms into equations (1) we ob- 
tain 

d d 
—m + {l,t) = - r c m + (l,t) + r r m~(l,t) - v + —m + (l,t) 

/OO 
m + (l',t)dl' (10a) 

d _ d 

—to (l,t) = + r c m + (l,t) - r r m (l,t)+v —to (l,t) 

at ol 

/OO 
[m + (l')+2m-(l')]dl', 

(10b) 

We note that in the special case when r c = r r =0, the 
equation for the length distribution of the growing micro- 
tubules is equivalent to the model of microtubule GTP- 
cap dynamics by Flyvbjerg et al. [20] with a vanishing 
diffusion constant. 

The number of parameters in these equations can be 
reduced by scaling the parameters and variables using 



natural units for length, time and microtubule number. 
As a unit of time, we take the mean time to catastrophe 
(l/r c ); the unit of length is the microtubule growth in 
that time (v + /r c ) and unit of microtubule number is the 
number of microtubules nucleated in this time (r n /r c ). 



In terms of the dimensionless length unit x 
define 
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(13) 
(14) 
(15) 

Note that the definition of v is inverted with respect to 
its natural conversion into dimensionless parameters, but 
this choice simplifies the notation in what follows. In the 
steady state, the equations (10) in dimensionless form 
become 

^-f + (x)=- f+(x)+rf-(x)~ S xf+(x) 
dx 
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f + (x')dx' (16a) 

J X 

f- (x) =-f + [x) + rf- (x) + sxf- (x) 

-s [f + (x') + 2f-(x')}dx\ (16b) 

J X 

with the boundary conditions 
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v dx 



/ + (0) = 1 
In the absence of severing (s 



lim f+(x) = ; lim f~(x) = 0. 

;— »oo x— too 

(17) 

0) these equations are 
solved by the exponential functions f + (x) = e~^ 1 ~ rv ) x 
and f~(x) = ve- ( - 1 ~ r ^ x [18]. Note that these solutions 
are only valid for rv < 1, when the average microtubule 
length is bounded. 

We proceed to analyze the steady state equations (16) 
in a number of steps. First, we determine global proper- 
ties of the length distributions. Subsequently, we derive 
an explicit expression for the microtubule length distri- 
butions in the special case in which rescues are absent 
(r = 0). We conclude with a numerical method using 
which the distribution can be calculated for arbitrary pa- 
rameter values. 



III. RESULTS 

A. Number of microtubules 

Multiplying equations (16) by x, integrating from to 
oo and subtracting the results yields 

/>oo />oo 

v / f + {x)dx = / r(x)dx. (18) 
Jo Jo 
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Also, integrating over (16a) yields 

/>oo 

- / + (0) = - / [f + (x) - rf-(x)] dx, (19) 
Jo 

which can be combined with (17) and (18) to give the 
following expression for the total number of microtubules. 



(f+( x ) + f-( x ))dx=^- (20) 
1 - rv 



Surprisingly, the total number of microtubules in the sys- 
tem does not depend on the rate of microtubule severing. 
Specifically, we note that severing can not prevent the 
diverging microtubule count for rv — > 1. 



B. In the absence of rescue events 

In the special case r = (no rescue events), equation 
(16a) for f + (x) decouples from (16b) and can be solved 
analytically. This solution can then be used to obtain 
an expression for f~(x) as well. We introduce the func- 
tions F ± (x) = — f ± (x') dx' , allowing us to rewrite 
equations (16) as 



a) /++/ 



^(^F+(x) + (l + sx)F+(x)) 







I_d_ 

v dx 2 



F~(x) = - 4~F + (x) +sx4-F~(x) 
dx dx 

+ sF + (x) + 2sF~{x) 



(21) 



(22) 



From equations (18) and (20) we find that F + (0) = — 1 
(because r — 0), and using the boundary condition 
d x F+(0) = /+(0) = 1 we solve (21) to obtain 
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Figure 1: Results in the absence of rescues (v = 1/2, r — 0). 
The choice v = 1/2 is an approximate value based on the in- 
vivo measurements in interphase plant cells reported by Vos 
et al. [21]. a) Length distributions for three different values of 
the severing rate: s — (black), s = 1 (dark gray) and s = 3 
(light gray) . The inset shows the same on a logarithmic scale, 
b) Average length as a function of the severing parameter s. 
The squares indicate the parameter values of the distributions 
in (a). The dashed line is the small- s approximation (A5). 
Length distributions have been calculated from (23) and (24); 
the average length from (29). 



F+( 3 



-exp 



-sx 



(23) 



Inserting this solution into equation (22) , it can be solved 
using the boundary conditions F~(0) = —v (from equa- 
tion (20)) and lim^^oo F~(0) = to give 



F (x) = — v exp 
1 



1 , 

-x sv' 

2 



erfc 



^^s(l + v)xexp 
1 + s(l + v)x 



(l + s(l + v)x)' 
2s(l + v) 



\ y/2s(l + V) 

where erfc(z) is the complementary error function 



erfc(y) 



-* 2 dt. 



(24) 



(25) 



The expressions for f + {x) and f~(x) follow by differenti- 
ation of (24). The resulting distribution for v =1/2 and 
various values of s is shown in figure 1. 



Looking at the derivative of the combined distribution 
f + (x) + f~(x), we find that for x = 



d.r 



(f + (x) + r(x))\ x=Q = -(1 + w)(l - s(l - 2v)). (26) 



The perhaps surprising implication is that if v < 1/2, 
sufficiently high severing rates can lead to a positive slope 
of the distribution at x — 0. In other words, in that case 
the length distribution no longer decreases monotonically 
with length. 

Finally, we compute the average microtubule length as 



f oo x(/+(x) + rGx))dx 

J Q °(f + (x) + f-(x))dx 



(27) 



1 + Wo 



(F+(x) + F-(x)) dx (28) 



z 



c~ l dt, 



l/z 



(29) 



with z = s(l +v). This confirms that the average length 
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decreases with increasing severing activity, a fact that is 
reflected in figure 1. 



C. Numerical evaluation of the length distribution 

For the more general case r > we use a numerical 
method to solve equations (16). In order to do so, we 
reparameterize the problem. Inspired by the results (18) 
and (20) we define 



p(x) = f+(x) + f-(x) 
q(x) = vf+(x)-f-(x). 



(30) 
(31) 



Hence, p(x) is a (dimensionless) measure of the total mi- 
crotubule density. In terms of these variables, the steady 
state equations (16) can be written as 



d t ^ 
ax 



-xp(x) 



p{x')dx' 



(32) 



dx 



p(x) = —(1 — rv)p(x) — (1 + r)q(x) 



s(l — v)xp(x) — sxq{x) + 

[q(x') + (1 - 2v)p{x')] Ax'. (33) 



Equation (32) can be formally solved to give 



q(x) = sv 



x p(x ) 



p(x")dx'' 



sv / [(2x-x')p(x')]dx'. 



dx' (34) 
(35) 



Inserting this expression into equation (33) we transform 
the problem into a single integral equation. We can sub- 
sequently remove the integrals by differentiating three 
times and obtain the linear fourth-order ODE 

p (4) 

(x) =(—1 + rv — s(l — v)x)p < - 3 ^ (x) 

+ s(-4 + v{5 + x(l + r + sx)))p (2) (x) 

+ Asv(l + r + 2sx)p {1) (x) + 12vs 2 p(x), (36) 

where p^ n '(x) stands for (d/dx) n p(x). To derive the 
boundary conditions for this problem we use (20), (18) 
and (17) to derive 



p{x) dx 



1 + v 



1-rv' 
q(x) dx = 0, 

q(0) = -p(0) + 1 + v. 



(37) 

(38) 
(39) 



By repeated application of these equalities and differenti- 
ation of equations (32) and (33) we obtain the boundary 



conditions 
p (1) (0) =(1 + v)rp(0) - (1 + v)(l + r) 



+ s(l + v) 



l-2v 
1 — rv 



(40) 



P (2) (0) = H"(l + v )0- - rv ) + 3s1 '] + (! + «) x 



(1 + r)(l — rv) — s 



3 — rv + 2rv i 



1 — rv 

p (3) (0) =r(l + v) [(1 - rv) 2 + sr(-3 + 7v)] p(0) 
- (1 + v)(l + r)(l - rv) 2 
+ s(l + v) [6 + 3r - 4v - hrv + 2rv 2 ] 



(41) 



1 + v 

1 — rv 



(3 - Av + 8v 2 ). 



(42) 



We note that these boundary conditions are of the form 



/ P(0) \ 
p^(0) 

\p {3) (o)J 



= P lP (0)+P 2 



(43) 



with a single undetermined parameter p(0). Because 
equation (36) is a homogeneous linear ODE, we can eval- 
uate it twice, using both Pi and P2 as boundary con- 
ditions. We generally find that both solutions diverge 
with opposing signs. The value of p(0) can therefore be 
determined from the constraint lim^^oo p(x) = 0. 

To this end, both solutions should be evaluated over 
a range that is as large as possible, whilst maintaining a 
very high numerical accuracy, because the final result is 
obtained by subtracting two diverging functions. In our 
numerical calculations (using Wolfram Mathcmatica 6.0), 
wc have evaluated the differential equation with sufficient 
precision to achieve an accuracy of 13 significant digits 
and limited the integration range to 10 decay lengths 
(.x = 10(1 — to) -1 ) or the point at which the first solu- 
tion exceeded the value 10 s , depending on which occurred 
first. 

Figure 2 shows the numerically computed distributions 
for v = 1/2, r = 1 and various values of s. It is interesting 
to note that the total distribution is no longer monoton- 
ically decreasing for s — 1 and s = 3. Figure 3 shows 
that this is solely due to the contribution from the grow- 
ing microtubules. We also note that the average length 
decreases rapidly for relatively small severing rate. As 
the severing rate increases, the average length converges 
to that of the system without rescue events (r = 0). In 
other words, if severing events occur very frequently, res- 
cue events are no longer significant, presumably because 
the microtubules become so short that they disappear be- 
fore they can be rescued. This statement is summarized 
by the condition 1/ry ^> l/v~, where I is the average 
microtubule length, or in dimensionless units rvx <^ 1. 
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Figure 2: Numerically computed length distributions and 
comparison with simulation data. Based on the in-vivo mea- 
surements in interphase plant cells reported by Vos et al. [21], 
approximate values of v = 1/2 and r = 1 have been used, 
a) Length distributions for three different values of the sev- 
ering rate: s — (black), s = 1 (dark gray) and s = 3 
(light gray). The inset shows the same on a logarithmic scale. 
Also shown for s = 1 is data obtained by a stochastic sim- 
ulation of microtubule severing using a Gillespie algorithm, 
matching the predicted distribution. Simulation parameters 
were [v + — 0.1 /im s" 1 , v~ — 0.2 /im s" 1 , r c — 0.01s~ , 
ry = 0.01s" 1 , r„ = 10s" 1 , r s = 0.001 pm _1 s -1 ]. Length 
data was distributed into 500 bins and sampled 1000 times 
at 50 second intervals after an initial equilibration period of 
50,000 seconds, b) Average length as a function of the sever- 
ing parameter s. The squares indicate the parameter values 
of the distributions in (a). The dashed line is the small- s ap- 
proximation (A5) The light gray curve, which converges for 
large s, is the result for r = (figure 1). Length distributions 
have been calculated from (36). 



IV. DISCUSSION 

We have constructed a model that describes the dy- 
namic instability of microtubules in combination with 
the severing of microtubules. This model takes the form 
of two coupled intcgro-diffcrential equations that are a 
function of three parameters: v, the ratio of polymeriza- 
tion and depolymerization speeds, r, the ratio of the res- 
cue and catastrophe rates, and the dimensionless severing 
rate s. We have analyzed the steady state solutions and 
their properties: notably the number of microtubules and 
their average length. For the special cases of no rescue 




Figure 3: Length distributions for growing and shrinking mi- 
crotubules separately (v = 1/2, r = 1, s = 1). The total 
distribution function p(x) = f + (x) + f~(x) has been com- 
puted numerically using (36); f + (x) and f~(x) have been 
computed using (35) and (31). 



events (r = 0) and small severing rates (see appendix A), 
we have presented analytical solutions. The general case 
has been addressed by transforming the coupled integro- 
differential equations into a single fourth order differen- 
tial equation that is solved numerically. The resulting 
microtubule length distributions show a number of inter- 
esting properties. 



1. Shorter, more compact length distributions 

As expected, an increase in the severing rate always 
leads to a decrease in the average length of the micro- 
tubules. In addition, because the rate of severing is 
proportional to the microtubule length, the number of 
very long microtubules is strongly reduced, leading to a 
distribution that is more compact (see appendix A for 
an explicit expression). Furthermore, in contrast to the 
dynamic instability model without severing, the length 
distributions are no longer guaranteed to decrease mono- 
tonically with increasing length. For some parameters, 
a 'bump' is observed in the length distribution, caused 
by the continuous creation of short (but not vanishingly 
small) microtubules through the severing process. 

2. Conservation of microtubule number 

We have also demonstrated that the total number of 
microtubules does not depend on the severing rate, even 
though the average length of each microtubule decreases. 
To understand this counterintuitive result, we consider a 
system without severing, in which the following steady- 
state relation holds: [population size] = [nucleation rate] 
x [average lifetime]. Suppose that a single microtubule 
is severed (consistent with an infinitesimally small sever- 
ing rate). The expected lifetime of the created segments 
is obviously shorter than that of the original segment. 
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However, the fact that the total number of microtubules 
is unaffected by the presence of severing implies that 
this decrease in lifetime is compensated exactly by the 
fact that every severing event 'nucleates' an extra micro- 
tubule. For this cancellation to occur, the sum of the 
lifetimes of the two fragments of a severed microtubule 
should equal the expected lifetime of the microtubule if 
severing had not occurred. 

That this is indeed the case can be shown using a sim- 
ple argument, based on the lifetime of a microtubule in 
the absence of severing. In this case the dynamics of 
the microtubule tip are independent of its length. This 
implies that the average time it takes for a growing mi- 
crotubule of length I to later return (in a shrinking state) 
to the same length is independent of the value of I, and 
therefore equal to the average microtubule lifetime T + 
(obtained by starting from l(to) = 0). The remaining 
contribution to the lifetime is given by the time it takes 
for the shrinking microtubule of length I to disappear. 
In the absence of rescue events, this would take a time 
l/v~. However, each rescue event switches the micro- 
tubule back to the growing state, extending the micro- 
tubule's expected lifetime by T + before it returns to the 
same position in the shrinking state. The expected num- 
ber of such rescue events is equal to r r l/v~. This argu- 
ment is illustrated in figure 4. Collecting the terms de- 
scribed above, we obtain the following expression for the 
expected microtubule lifetime T(l, a) of a microtubule of 
length I and growth state a £ { — ,+}: 

T(l,a)=T+5 a ,+ + — (l + r r T+). (44) 
v 

The same relation, including an expression for T + , has 
been derived analytically by and Bicout [22] and an equiv- 
alent expression for microtubules that grow in discrete 
units was produced by Rubin [23]. 

Each severing event conserves the total length I of the 
severed microtubule and the state a of the existing plus 
end, and the newly created plus end immediately under- 
goes a catastrophe (<r new = — )• In terms of life times we 
find that T(l,a) = T(l', a) + T(l — I', — ), so that a single 
severing event indeed preserves the total lifetime. By in- 
duction it follows that any number of severing events will 
leave the expected number of microtubules in the system 
intact. 

We stress that this argument holds regardless of the fre- 
quency and location of severing events. This implies that 
any distribution of severing events - provided they lead 
to a catastrophe of the trailing end - conserves the total 
number of microtubules in steady state. Specifically, this 
also applies to severing at positions where microtubules 
of different orientations cross [14] , a process that may be 
relevant to the formation of the cortical array in plant 
cells. 

The conservation of microtubule number in the pres- 
ence of severing is in apparent contradiction with the ex- 
perimentally reported increase in microtubule numbers 
[11]. However, this can be readily explained by the fact 
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time 

Figure 4: Schematic depiction of the life history of a shrinking 
microtubule with an initial length 1(0). The return times U 
have an expectation value T + . 



that our results are based on the assumption of constant 
parameter values. In a living cell, the decrease in aver- 
age microtubule length that is the result of severing will 
lead to an increased availability of free tubulin dimcrs. 
In turn, this is likely to increase the polymerization rate 
(growth speed) and nucleation rate, which would indeed 
cause an increase in the number of microtubules. 



Acknowledgments 

ST thanks Frank Jiilicher and Benjamin Lindner for 
very helpful discussions and is grateful for the hospitality 
that was provided by the Max Planck Institute for the 
Physics of Complex Systems. The authors thank Nils 
Becker for critical comments on the manuscript. ST was 
supported by a grant from the NWO "Computational 
Life Sciences" program (contract CLS 635.100.003) and a 
visiting stipend from the Max Planck Society. This work 
is part of the research program of the "Stichting voor Fun- 
damcntccl Ondcrzock dor Materie (FOM)", which is fi- 
nancially supported by the "Nederlandsc organisatie voor 
Wetenschappelijk Onderzoek (NWO)". 



Appendix A: Small severing rates 

In the absence of severing, equations (16) are solved 
by f + (x) = c - (1 - rv ^ x and f~(x) = we^ 1 "™^ . To inves- 
tigate the changes that occur when a small severing rate 
is included, we perturb these solutions as follows 

/+(*) = c- {1 - rv)x {l + sf+(x)) + 0(s 2 ) (Al) 
f-(x) = ve- {1 ~ rv)x (l + sf-(x)) + 0(s 2 ). (A2) 

Inserting these expressions into equations (16) and drop- 
ping all higher order terms gives a set of equations that 
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can be solved to yield 



Finally, we determine the length variation 



f+( x ) =e -(i-«0* h+ a 

+ 0{s 2 ) 
f-(x) =ve- {1 - rv)x (l + s 



1 + rv 1 + rv 



{1-rv) 2 2(1 -rv) 

(A3) 

1 + v —v + rv + 2rv 2 



1 + rv 2 2 
2(1 - rv)^ 



(1 — rv) 2 (1 — rv) 2 

+ 0(s 2 ) (A4) 



We note that both solutions will become negative for 
large values of x. However, even though this is decidedly 
unphysical, the effect on measurable parameters such as 
the average microtubule length is small (for small s), be- 
cause of the rapid decay of |/ + (a;)| and |/~(a;)|. The 
properties of the resulting distributions, such as the av- 
erage length, will therefore still remain valid, subject to 
the bounds on s computed below. 

The first-order distributions in s satisfy the total mi- 
crotubule number constraint J Q [f + (x) + f~(x)]dx = 
(1 + v)/(l — rv) + 0(s 2 ), consistent with the general re- 
sult (20). For the average length we obtain (to first order 
in s) 



1 



1 



1 + v 

'(1 - rv) 4 



0(s 2 



(A5) 



For r = 0, this is consistent with the exact result (29). 
The resulting (linear) predictions are indicated in figures 
1 and 2. 



u% =< (x-xY >= - 



(1 — TVY 



(1 — rv) 5 



and, from that, the coefficient of variation (o~ x /x) 



(A6) 



1 _|_ T*7j2 

a x /x = l-k - +rv +0(s 2 ). (A7) 
(1 — rv) 6 

This number provides a measure for the relative width 
of the distribution. The results (A5) and (A7) confirm 
that severing decreases both the weighted average and 
the relative width of the length distribution. 

The results above have been obtained under the as- 
sumption that s is very small. To make an a priori esti- 
mate for the validity range of s, we determine the relative 
importance of the terms on the right-hand side of equa- 
tions (16a) and (16b). Using the results in the absence of 
severing as a benchmark, the terms not involving s give 
contributions of the order (1 — rv)e~( 1 ~ rv ^ x . Comparing 
the integral term (evaluated for the s = situation) with 
this term gives s <C (1 — rv) 2 . The term that is pro- 
portional to sx will dominate the other terms for large 
x, but this does not significantly affect the results if it 
only occurs for lengths that are much longer than the 
average length. Evaluating the terms at x = n/(l — rv), 
where n is the number of average lengths, we obtain the 
constraint s <C (1 — rv) /n. Because n is of the order 
1 , we simply state that the approximation is accurate for 
s < (1 - rv) 2 . 
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